<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.9.1"/>
<title>hybridMPIC: particles.h Source File</title>
<link href="../../tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="../../jquery.js"></script>
<script type="text/javascript" src="../../dynsections.js"></script>
<link href="../../search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="../../search/searchdata.js"></script>
<script type="text/javascript" src="../../search/search.js"></script>
<script type="text/javascript">
  $(document).ready(function() { init_search(); });
</script>
<link href="../../doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  <td id="projectlogo"><img alt="Logo" src="../../logo.png"/></td>
  <td style="padding-left: 0.5em;">
   <div id="projectname">hybridMPIC
   </div>
   <div id="projectbrief">A fully relativisitic 2d Particle-In-Cell(PIC) code, with explicit solver, MPI&amp;OpenMP support</div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.9.1 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "../../search",false,'Search');
</script>
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="../../index.html"><span>Main&#160;Page</span></a></li>
      <li><a href="../../pages.html"><span>Related&#160;Pages</span></a></li>
      <li><a href="../../annotated.html"><span>Data&#160;Structures</span></a></li>
      <li class="current"><a href="../../files.html"><span>Files</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
        <span class="left">
          <img id="MSearchSelect" src="../../search/mag_sel.png"
               onmouseover="return searchBox.OnSearchSelectShow()"
               onmouseout="return searchBox.OnSearchSelectHide()"
               alt=""/>
          <input type="text" id="MSearchField" value="Search" accesskey="S"
               onfocus="searchBox.OnSearchFieldFocus(true)" 
               onblur="searchBox.OnSearchFieldFocus(false)" 
               onkeyup="searchBox.OnSearchFieldChange(event)"/>
          </span><span class="right">
            <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="../../search/close.png" alt=""/></a>
          </span>
        </div>
      </li>
    </ul>
  </div>
  <div id="navrow2" class="tabs2">
    <ul class="tablist">
      <li><a href="../../files.html"><span>File&#160;List</span></a></li>
    </ul>
  </div>
</div><!-- top -->
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
     onmouseover="return searchBox.OnSearchSelectShow()"
     onmouseout="return searchBox.OnSearchSelectHide()"
     onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>

<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0" 
        name="MSearchResults" id="MSearchResults">
</iframe>
</div>

<div class="header">
  <div class="headertitle">
<div class="title">particles.h</div>  </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="preprocessor">#ifndef __particles</span></div>
<div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="preprocessor">#define __particles</span></div>
<div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="preprocessor">#include &lt;iostream&gt;</span></div>
<div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="preprocessor">#include &lt;vector&gt;</span></div>
<div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="preprocessor">#include &quot;type.h&quot;</span></div>
<div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="preprocessor">#include &quot;physicalconstants.h&quot;</span></div>
<div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="preprocessor">#include &quot;stdio.h&quot;</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html">    9</a></span>&#160;<span class="keyword">class </span><a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>{</div>
<div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;    <span class="keyword">public</span>:</div>
<div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;        <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>(<a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; pos, <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a>&amp; vel, <span class="keywordtype">long</span>&amp; wt,\</div>
<div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;                <span class="keywordtype">double</span>&amp; ene, <span class="keywordtype">double</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a>):<a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a>(pos), <a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a>(pos), <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">\</a></div>
<div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">                velocity</a>(vel), <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a>(vel), <a class="code" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">weight</a>(wt), <a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a>(ene), <a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a>(charge){</div>
<div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;                <span class="keywordtype">id</span> = <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a>;</div>
<div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;                <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a> ++;</div>
<div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;                <a class="code" href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">total_num</a> ++;</div>
<div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;                lorentz_gamma = 1.0 / sqrt(1.0 - (<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> * <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>) / (c_light_speed\</div>
<div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;                            * c_light_speed));</div>
<div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;                <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> zero(0.0, 0.0, 0.0);</div>
<div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;                <a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a> = zero;</div>
<div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;                <a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a> = zero;</div>
<div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;        }</div>
<div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;        <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>(<span class="keyword">const</span> <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>&amp; other){</div>
<div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a190602549e77f92d18694bfac93774de">species</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a190602549e77f92d18694bfac93774de">species</a>;</div>
<div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a>;</div>
<div class="line"><a name="l00026"></a><span class="lineno">   26</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>;</div>
<div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a>;</div>
<div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a>;</div>
<div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;            this-&gt;charge = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a>;</div>
<div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;            this-&gt;<span class="keywordtype">id</span> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#aeba79577ec36efdd5d55f0b483663eeb">id</a>;</div>
<div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;            <span class="comment">//cout&lt;&lt;&quot;in ptcls copy func &quot;&lt;&lt;other.id&lt;&lt;endl;</span></div>
<div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">weight</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">weight</a>;</div>
<div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a>;</div>
<div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a>;</div>
<div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a>;</div>
<div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;            this-&gt;<a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a> = other.<a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a>;</div>
<div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;            this-&gt;lorentz_gamma = other.lorentz_gamma;</div>
<div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">total_num</a> ++;</div>
<div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;        }</div>
<div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;        <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>(){</div>
<div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;            <span class="keywordtype">id</span> = <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a>;</div>
<div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a> ++;</div>
<div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">total_num</a> ++;</div>
<div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;            <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> zero(0.0, 0.0, 0.0);</div>
<div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> = zero;</div>
<div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a> = zero;</div>
<div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;            lorentz_gamma = 1.0;</div>
<div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a> = 0.0;</div>
<div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a> = zero;</div>
<div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;            <a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a> = zero;</div>
<div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;        }</div>
<div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;</div>
<div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;        <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>(<span class="keyword">const</span> <a class="code" href="../../d0/de9/structParticles.html">Particles</a>&amp; myptcl);</div>
<div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;    <span class="keyword">public</span>:</div>
<div class="line"><a name="l00056"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">   56</a></span>&#160;        <span class="keyword">static</span> <span class="keywordtype">long</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a>;</div>
<div class="line"><a name="l00058"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">   58</a></span>&#160;        <span class="keyword">static</span> <span class="keywordtype">long</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">total_num</a>;</div>
<div class="line"><a name="l00060"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#aeba79577ec36efdd5d55f0b483663eeb">   60</a></span>&#160;        <span class="keywordtype">long</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#aeba79577ec36efdd5d55f0b483663eeb">id</a>;</div>
<div class="line"><a name="l00062"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a190602549e77f92d18694bfac93774de">   62</a></span>&#160;        <span class="keywordtype">int</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a190602549e77f92d18694bfac93774de">species</a>;</div>
<div class="line"><a name="l00064"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">   64</a></span>&#160;        <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a>;</div>
<div class="line"><a name="l00066"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">   66</a></span>&#160;        <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a>;</div>
<div class="line"><a name="l00068"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">   68</a></span>&#160;        <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>;</div>
<div class="line"><a name="l00070"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">   70</a></span>&#160;        <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a>;</div>
<div class="line"><a name="l00072"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">   72</a></span>&#160;        <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a>;</div>
<div class="line"><a name="l00074"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">   74</a></span>&#160;        <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> <a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a>;</div>
<div class="line"><a name="l00076"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">   76</a></span>&#160;        <span class="keywordtype">long</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">weight</a>;</div>
<div class="line"><a name="l00078"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">   78</a></span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a>;</div>
<div class="line"><a name="l00080"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">   80</a></span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a>;</div>
<div class="line"><a name="l00082"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">   82</a></span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a>;</div>
<div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;        <span class="keywordtype">double</span> lorentz_gamma;</div>
<div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;        <span class="keywordtype">double</span> od_factor;</div>
<div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;        <span class="keywordtype">double</span> rod_factor;</div>
<div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;        <span class="keywordtype">double</span> qed_factor;</div>
<div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;        ~<a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>(){</div>
<div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;            total_num --;</div>
<div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;        }</div>
<div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;        <span class="keyword">friend</span> ostream&amp; operator&lt;&lt;(ostream&amp; out, <span class="keyword">const</span> <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>&amp; myptcl)</div>
<div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;        {</div>
<div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;            out&lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#aeba79577ec36efdd5d55f0b483663eeb">id</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>&lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>&lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>&lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>\</div>
<div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;                &lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">weight</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>&lt;&lt;myptcl.<a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a>&lt;&lt;<span class="stringliteral">&quot;\t&quot;</span>&lt;&lt;myptcl.lorentz_gamma;</div>
<div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;            <span class="keywordflow">return</span> out;</div>
<div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;        }</div>
<div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#af863f11743d992891e0e79de66e28f65">advance</a>(<span class="keywordtype">double</span>&amp; delta_t);</div>
<div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a522b307a9a3fefda5c1cf489ca77fa9b">boris</a>(<span class="keywordtype">double</span>&amp; delta_t);</div>
<div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a18c3f527195910edacebf07d1deee7c3">lorentz</a>(<span class="keywordtype">double</span>&amp; delta_t);</div>
<div class="line"><a name="l00103"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a77ff6b30069c0a76772206cebceb009f">  103</a></span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a77ff6b30069c0a76772206cebceb009f">set_field</a>(<span class="keyword">const</span> <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a>&amp; oefield, <span class="keyword">const</span> <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a>&amp; obfield){</div>
<div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;            efield = oefield;</div>
<div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;            bfield = obfield;</div>
<div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;        }</div>
<div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;</div>
<div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;};</div>
<div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;</div>
<div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;</div>
<div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;</div>
<div class="line"><a name="l00112"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#af863f11743d992891e0e79de66e28f65">  112</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#af863f11743d992891e0e79de66e28f65">__Ptcls::advance</a>(<span class="keywordtype">double</span>&amp; delta_t){</div>
<div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a522b307a9a3fefda5c1cf489ca77fa9b">boris</a>(delta_t);</div>
<div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a18c3f527195910edacebf07d1deee7c3">lorentz</a>(delta_t);</div>
<div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;}</div>
<div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;</div>
<div class="line"><a name="l00117"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a18c3f527195910edacebf07d1deee7c3">  117</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a18c3f527195910edacebf07d1deee7c3">__Ptcls::lorentz</a>(<span class="keywordtype">double</span>&amp; delta_t){</div>
<div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">old_pos</a> = <a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a>;</div>
<div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;    <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a> vel(<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> + <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a>);</div>
<div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a> = <a class="code" href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">position</a> + vel * delta_t * 0.5;</div>
<div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;    <span class="comment">//cout&lt;&lt;id&lt;&lt;&quot;\t&quot;&lt;&lt;species&lt;&lt;&quot; pos &quot;&lt;&lt;position.member[0]&lt;&lt;&quot; old &quot;&lt;&lt;old_pos.member[0]&lt;&lt;endl;</span></div>
<div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;    <span class="comment">//getchar();</span></div>
<div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">old_vel</a> = <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>;</div>
<div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;}</div>
<div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;</div>
<div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;</div>
<div class="line"><a name="l00127"></a><span class="lineno"><a class="line" href="../../d9/ddf/class____Ptcls.html#a522b307a9a3fefda5c1cf489ca77fa9b">  127</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="../../d9/ddf/class____Ptcls.html#a522b307a9a3fefda5c1cf489ca77fa9b">__Ptcls::boris</a>(<span class="keywordtype">double</span>&amp; delta_t){</div>
<div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;    <span class="comment">//cout&lt;&lt;&quot; pid &quot;&lt;&lt;id&lt;&lt;&quot; ps &quot;&lt;&lt;species&lt;&lt;endl;</span></div>
<div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a> == 0)</div>
<div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;    {</div>
<div class="line"><a name="l00131"></a><span class="lineno">  131</span>&#160;        <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;    }</div>
<div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> zero(0, 0, 0);</div>
<div class="line"><a name="l00134"></a><span class="lineno">  134</span>&#160;    lorentz_gamma = 1.0 / sqrt(1.0 -(<a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> * <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>) / c_light_speed / c_light_speed);</div>
<div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160;    <span class="keywordtype">double</span> qtm2 = delta_t * <a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a> * 0.5 / <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a>;</div>
<div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> uinit = lorentz_gamma * <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a>;</div>
<div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> uminus = uinit + qtm2 * <a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a>;</div>
<div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;    </div>
<div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;    lorentz_gamma = sqrt(1.0 + (uminus * uminus) / c_light_speed / c_light_speed);</div>
<div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> tvec = qtm2 * <a class="code" href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">bfield</a> / lorentz_gamma;</div>
<div class="line"><a name="l00141"></a><span class="lineno">  141</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> uslash = uminus + Cross(uminus, tvec);</div>
<div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> svec = 2.0 * tvec / (1.0 + tvec * tvec);</div>
<div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> uplus = uminus + Cross(uslash, svec);</div>
<div class="line"><a name="l00144"></a><span class="lineno">  144</span>&#160;    </div>
<div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> ufinal = uplus + qtm2 * <a class="code" href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">efield</a>;</div>
<div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;    lorentz_gamma = sqrt(1.0 + (ufinal * ufinal) / c_light_speed / c_light_speed);</div>
<div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;    velocity = ufinal / lorentz_gamma;</div>
<div class="line"><a name="l00148"></a><span class="lineno">  148</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a> = (lorentz_gamma - 1.0) * <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a> * c_light_speed * c_light_speed;</div>
<div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;    <span class="keywordflow">if</span>(<a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a> &lt; 0)</div>
<div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;    {</div>
<div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;        cout&lt;&lt;<span class="stringliteral">&quot; negative energy !&quot;</span>&lt;&lt;endl;</div>
<div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;    }</div>
<div class="line"><a name="l00153"></a><span class="lineno">  153</span>&#160;    </div>
<div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;    <span class="keywordflow">if</span>(Mod(velocity)&gt;c_light_speed)</div>
<div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;    {</div>
<div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;        cout&lt;&lt;<span class="stringliteral">&quot; velocity error &quot;</span>&lt;&lt;efield&lt;&lt;endl;</div>
<div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;    }</div>
<div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;</div>
<div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;}</div>
<div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;</div>
<div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;</div>
<div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;</div>
<div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;<span class="comment">/***</span></div>
<div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;<span class="comment"> * how to cycle this ptclslist?</span></div>
<div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;<span class="comment"> *      1st declare an list and</span></div>
<div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;<span class="comment"> * typedef vector&lt;vector&lt;__Electron&gt; &gt; mylist;</span></div>
<div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;<span class="comment"> *      2nd declare an iterator</span></div>
<div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;<span class="comment"> * vector&lt;__Electron&gt;::iterator myiter;</span></div>
<div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;<span class="comment"> *      3rd use this method to cycle and</span></div>
<div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;<span class="comment"> * for (myiter = myelist.begin(); myiter != myelist.end(); myiter ++)</span></div>
<div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;<span class="comment"> * {</span></div>
<div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;<span class="comment"> *      cout&lt;&lt;&quot;id of my ptcls&quot;&lt;&lt;myiter-&gt;id&lt;&lt;endl;</span></div>
<div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;<span class="comment"> *      }</span></div>
<div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;<span class="comment"> * I think i should write a iterator, so that it will be more flexiable</span></div>
<div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;<span class="comment"> *      ***/</span></div>
<div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;</div>
<div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;</div>
<div class="line"><a name="l00182"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html">  182</a></span>&#160;<span class="keyword">class </span><a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>{</div>
<div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;    <span class="keyword">public</span>:</div>
<div class="line"><a name="l00185"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">  185</a></span>&#160;        vector&lt;__Ptcls&gt; <a class="code" href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">member</a>;</div>
<div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;        <span class="comment">// use this func you can call by obj.iter() to get a iter</span></div>
<div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        <span class="comment">// then use it to find any particle you want</span></div>
<div class="line"><a name="l00189"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a98e6733c424414f7cc596e64486e43b2">  189</a></span>&#160;<span class="comment"></span>        <span class="keyword">typename</span> vector&lt;__Ptcls&gt;::iterator <a class="code" href="../../d6/d19/class____Ptclslist.html#a98e6733c424414f7cc596e64486e43b2">iter</a>;</div>
<div class="line"><a name="l00191"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a639c350be512a5b80043429e5e873a75">  191</a></span>&#160;        <span class="keywordtype">int</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a639c350be512a5b80043429e5e873a75">species</a>;</div>
<div class="line"><a name="l00193"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#ac6fc791bb5b3b456b7a909c175f43340">  193</a></span>&#160;        <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a> <a class="code" href="../../d6/d19/class____Ptclslist.html#ac6fc791bb5b3b456b7a909c175f43340">mylength</a>;</div>
<div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;        <span class="comment">// 0 e- 1 i+ 2 e+ 4 gamma</span></div>
<div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;        <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>(){}</div>
<div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;        <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; mylength, <span class="keywordtype">int</span> sp):mylength(mylength), species(sp){}</div>
<div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;        <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>(<span class="keyword">const</span> <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>&amp; other){</div>
<div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;            this-&gt;member = other.<a class="code" href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">member</a>;</div>
<div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;            this-&gt;species = other.<a class="code" href="../../d6/d19/class____Ptclslist.html#a639c350be512a5b80043429e5e873a75">species</a>;</div>
<div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;            this-&gt;mylength = other.<a class="code" href="../../d6/d19/class____Ptclslist.html#ac6fc791bb5b3b456b7a909c175f43340">mylength</a>;</div>
<div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;        }</div>
<div class="line"><a name="l00203"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a682e2a2543c54e627c6fc3fd42843878">  203</a></span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a682e2a2543c54e627c6fc3fd42843878">squeeze</a>()</div>
<div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;        {</div>
<div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;            <span class="keywordflow">if</span> (member.capacity() &gt; 2 * member.size())</div>
<div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;            {</div>
<div class="line"><a name="l00207"></a><span class="lineno">  207</span>&#160;                vector&lt;__Ptcls&gt;(<a class="code" href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">member</a>).swap(member);</div>
<div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;            }</div>
<div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;        }</div>
<div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a67f83b452b3c0970072c6b8454fe8835">merge</a>(<span class="keyword">const</span> <span class="keywordtype">double</span>&amp; delta_e, <span class="keyword">const</span> <span class="keywordtype">double</span>&amp; delta_v);</div>
<div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;        <span class="comment">// use this push func, different species ptcls cannot be</span></div>
<div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;        <span class="comment">// stored in same vector, since the are different __Ptcls(ty</span></div>
<div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;        <span class="comment">// pe)</span></div>
<div class="line"><a name="l00216"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a7969ef614256281a5103a4478e009e46">  216</a></span>&#160;<span class="comment"></span>        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a7969ef614256281a5103a4478e009e46">push</a>(<span class="keyword">const</span> <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>&amp; other)</div>
<div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;        {</div>
<div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;            member.push_back(other);</div>
<div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;            <span class="comment">//cout&lt;&lt;&quot;in pushing -----------------------------------&quot;&lt;&lt;endl;</span></div>
<div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;            <span class="comment">//cout&lt;&lt;other.lorentz_gamma&lt;&lt;endl;</span></div>
<div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;            <span class="comment">//cout&lt;&lt;member.size()&lt;&lt;endl;</span></div>
<div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;            <span class="comment">//cout&lt;&lt;&quot;pushin over-----------------------------&quot;&lt;&lt;endl;</span></div>
<div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;            <span class="comment">//getchar();</span></div>
<div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160;        }</div>
<div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;        <span class="comment">// delete particle by offset, i.e. index, this is quicker</span></div>
<div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;        <span class="comment">// than other method</span></div>
<div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;        <span class="comment">// how to call?</span></div>
<div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;        <span class="comment">// just pop(n) n is the offset</span></div>
<div class="line"><a name="l00230"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a51c29b9954cc6f3fe9c51fb919d38a6d">  230</a></span>&#160;<span class="comment"></span>        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c29b9954cc6f3fe9c51fb919d38a6d">pop</a>(<span class="keywordtype">int</span> offset)</div>
<div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;        {</div>
<div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;            <span class="comment">//cout&lt;&lt;&quot;int plist offset&quot;&lt;&lt;endl;</span></div>
<div class="line"><a name="l00233"></a><span class="lineno">  233</span>&#160;            <span class="keywordflow">if</span> (offset &lt; member.size())</div>
<div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;            {</div>
<div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;                iter = member.begin();</div>
<div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;                std::advance(iter, offset);</div>
<div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;                member.erase(iter);</div>
<div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;            }</div>
<div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;        }</div>
<div class="line"><a name="l00241"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a89bd5130d9b94c556e1f23ccf1a98485">  241</a></span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a89bd5130d9b94c556e1f23ccf1a98485">clear</a>()</div>
<div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;        {</div>
<div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;            clear_vec(member);</div>
<div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;            <span class="comment">//member.clear();</span></div>
<div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;        }</div>
<div class="line"><a name="l00247"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a7feae0ba57fcfc32edaffa105cb1d09c">  247</a></span>&#160;        <span class="keywordtype">int</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a7feae0ba57fcfc32edaffa105cb1d09c">size</a>(){</div>
<div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;            <span class="keywordflow">return</span> member.size();</div>
<div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;        }</div>
<div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;        <span class="comment">// use this func, you can get the begin of member, no longer</span></div>
<div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;        <span class="comment">// need to care about what kind of its type</span></div>
<div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;        <span class="keyword">typename</span> vector&lt;__Ptcls&gt;::iterator begin(){</div>
<div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;            <span class="keywordflow">return</span> member.begin();</div>
<div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;        }</div>
<div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;        <span class="keyword">typename</span> vector&lt;__Ptcls&gt;::iterator end(){</div>
<div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;            <span class="keywordflow">return</span> member.end();</div>
<div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;        }</div>
<div class="line"><a name="l00259"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#ab01ad08b2b0eba942ccb4e7af6b05ea6">  259</a></span>&#160;        <a class="code" href="../../d9/ddf/class____Ptcls.html">__Ptcls</a>&amp; <a class="code" href="../../d6/d19/class____Ptclslist.html#ab01ad08b2b0eba942ccb4e7af6b05ea6">get_ptcl</a>(<span class="keywordtype">int</span> offset)</div>
<div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;        {</div>
<div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;            std::advance(iter, offset - 1);</div>
<div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;            <span class="keywordflow">return</span> *(<a class="code" href="../../d6/d19/class____Ptclslist.html#a98e6733c424414f7cc596e64486e43b2">iter</a>);</div>
<div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;        }</div>
<div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#abd6d139f042a71a595aa7f4aa2a21714">get_sp_ek</a>(<span class="keyword">const</span> <span class="keywordtype">int</span>&amp; sp);</div>
<div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;        ~<a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>(){</div>
<div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;            clear_vec(member);</div>
<div class="line"><a name="l00268"></a><span class="lineno">  268</span>&#160;            <span class="comment">//member.clear();</span></div>
<div class="line"><a name="l00269"></a><span class="lineno">  269</span>&#160;            this-&gt;<a class="code" href="../../d6/d19/class____Ptclslist.html#a682e2a2543c54e627c6fc3fd42843878">squeeze</a>();</div>
<div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;        }</div>
<div class="line"><a name="l00272"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#abf8f9068e29bee8f2cdef145ba7e32ff">  272</a></span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#abf8f9068e29bee8f2cdef145ba7e32ff">advance</a>(<span class="keywordtype">double</span> delta_t, vector&lt;__Ptcls&gt;&amp; lists, \</div>
<div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;                <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; xymin, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; xymax,\</div>
<div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;                <span class="keyword">const</span> <span class="keywordtype">int</span>&amp; myrank){</div>
<div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;            <span class="keywordflow">if</span>(member.size() == 0)</div>
<div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;            {</div>
<div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;                <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;            }</div>
<div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;            <span class="keywordtype">int</span> i = member.size() - 1;</div>
<div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;            <span class="keywordflow">while</span>(i &gt;= 0){</div>
<div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;                member[i].advance(delta_t);</div>
<div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;                <span class="keywordflow">if</span>((member[i].position.member[0] &lt; xymin.member[0] || \</div>
<div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;                            member[i].position.member[0] &gt; xymax.member[0] || \</div>
<div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;                            member[i].position.member[1] &lt; xymin.member[1] || \</div>
<div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;                            member[i].position.member[1] &gt; xymax.member[1]))</div>
<div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;                {</div>
<div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;                    lists.push_back(member[i]);</div>
<div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;                    <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c29b9954cc6f3fe9c51fb919d38a6d">pop</a>(i);</div>
<div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;                }</div>
<div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;                i --;</div>
<div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;            }</div>
<div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;        }</div>
<div class="line"><a name="l00294"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#aad813503cf966e45b694006e71a022e1">  294</a></span>&#160;        <span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#aad813503cf966e45b694006e71a022e1">showme</a>(){</div>
<div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;            <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i = 0; i&lt; member.size();i ++){</div>
<div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;                cout&lt;&lt;member[i]&lt;&lt;endl;</div>
<div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;            }</div>
<div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;        }</div>
<div class="line"><a name="l00300"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a86188e07e2f77f26cb58aaab4fa92743">  300</a></span>&#160;        <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a> <a class="code" href="../../d6/d19/class____Ptclslist.html#a86188e07e2f77f26cb58aaab4fa92743">operator = </a>(<span class="keyword">const</span> <a class="code" href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a>&amp; other){</div>
<div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;            this-&gt;member = other.<a class="code" href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">member</a>;</div>
<div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;            this-&gt;species = other.<a class="code" href="../../d6/d19/class____Ptclslist.html#a639c350be512a5b80043429e5e873a75">species</a>;</div>
<div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;        }</div>
<div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">weightfunc</a>(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; mycoords, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; other, \</div>
<div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;                <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; cell_length);</div>
<div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;        <span class="keywordtype">double</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">weightfunc</a>(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; mycoords, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; other);</div>
<div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;        <span class="keywordtype">double</span> weight_to_coords_rho(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; coords, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp;\</div>
<div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;        cell_length);</div>
<div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;        <span class="keywordtype">double</span> weight_to_coords_j(<span class="keyword">const</span> <span class="keywordtype">int</span>&amp; sp, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; coord, \</div>
<div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;                <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; cell_length);</div>
<div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;        <span class="comment">// 这里的sp指的是0---jx， 1--jy，2--jz</span></div>
<div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;};</div>
<div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;</div>
<div class="line"><a name="l00315"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">  315</a></span>&#160;<span class="keywordtype">double</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">__Ptclslist::weightfunc</a>(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; mycoords, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; other, \</div>
<div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;        <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; cell_length){</div>
<div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160;    <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a> temp(0, 0);</div>
<div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;    temp = (mycoords - other).Abs();</div>
<div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;    <span class="comment">//cout&lt;&lt;&quot;temp &quot;&lt;&lt;temp&lt;&lt;endl;</span></div>
<div class="line"><a name="l00320"></a><span class="lineno">  320</span>&#160;       <span class="keywordflow">if</span>(cell_length &gt; temp){</div>
<div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;           <span class="keywordflow">return</span> ((cell_length.member[0] - temp.member[0])\</div>
<div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160;                   * (cell_length.member[1] - temp.member[1]))\</div>
<div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;               / (cell_length.member[0] * cell_length.member[1]);</div>
<div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;       }</div>
<div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;    <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;        <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;    }</div>
<div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;}</div>
<div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;</div>
<div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;</div>
<div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;</div>
<div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;<span class="keywordtype">double</span> __Ptclslist::weight_to_coords_rho(<span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; coords, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp;\</div>
<div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;        cell_length){</div>
<div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;    <span class="keywordtype">double</span> rho = 0;</div>
<div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;    <span class="keywordflow">if</span>(member.size() != 0){</div>
<div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;        <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i = 0; i &lt; member.size(); i ++){</div>
<div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;            rho += member[i].weight * member[i].charge * <a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">\</a></div>
<div class="line"><a name="l00338"></a><span class="lineno">  338</span>&#160;<a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">                   weightfunc</a>(member[i].position, coords, cell_length);</div>
<div class="line"><a name="l00339"></a><span class="lineno">  339</span>&#160;</div>
<div class="line"><a name="l00340"></a><span class="lineno">  340</span>&#160;        }</div>
<div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;    }</div>
<div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;    <span class="keywordflow">return</span> rho;</div>
<div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;}</div>
<div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;</div>
<div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;<span class="keywordtype">double</span> __Ptclslist::weight_to_coords_j(<span class="keyword">const</span> <span class="keywordtype">int</span>&amp; sp, <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; coord, \</div>
<div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;        <span class="keyword">const</span> <a class="code" href="../../d6/d9c/class____Vect2.html">__Vect2&lt;double&gt;</a>&amp; cell_length){</div>
<div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;    <span class="keywordtype">double</span> rets = 0;</div>
<div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;    <span class="keywordflow">if</span>(member.size()!=0){</div>
<div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;        <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i = 0; i &lt; member.size(); i ++){</div>
<div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;            rets += member[i].weight * member[i].velocity.member[sp] *<a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">\</a></div>
<div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;<a class="code" href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">                    weightfunc</a>(member[i].position, coord, cell_length) * member[i].charge;</div>
<div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;        }</div>
<div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;    }</div>
<div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;    <span class="keywordflow">return</span> rets;</div>
<div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;}</div>
<div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;</div>
<div class="line"><a name="l00357"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#abd6d139f042a71a595aa7f4aa2a21714">  357</a></span>&#160;<span class="keywordtype">double</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#abd6d139f042a71a595aa7f4aa2a21714">__Ptclslist::get_sp_ek</a>(<span class="keyword">const</span> <span class="keywordtype">int</span>&amp; sp)</div>
<div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;{</div>
<div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;    <span class="keywordflow">if</span>(sp != species)</div>
<div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;    {</div>
<div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;        <span class="keywordflow">return</span> 0;</div>
<div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;    }</div>
<div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;    <span class="keywordtype">double</span> retek = 0;</div>
<div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;    <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i = 0; i &lt; member.size(); i ++)</div>
<div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;    {</div>
<div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;        retek += member[i].weight * member[i].energy;</div>
<div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;    }</div>
<div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;    <span class="keywordflow">return</span> retek;</div>
<div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;}</div>
<div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;</div>
<div class="line"><a name="l00371"></a><span class="lineno"><a class="line" href="../../d6/d19/class____Ptclslist.html#a67f83b452b3c0970072c6b8454fe8835">  371</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="../../d6/d19/class____Ptclslist.html#a67f83b452b3c0970072c6b8454fe8835">__Ptclslist::merge</a>(<span class="keyword">const</span> <span class="keywordtype">double</span>&amp; delta_e, <span class="keyword">const</span> <span class="keywordtype">double</span>&amp; delta_v)</div>
<div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;{</div>
<div class="line"><a name="l00373"></a><span class="lineno">  373</span>&#160;    <span class="keywordtype">int</span> nid = 0;</div>
<div class="line"><a name="l00374"></a><span class="lineno">  374</span>&#160;    <a class="code" href="../../da/d1a/class____Vect3.html">__Vect3&lt;double&gt;</a> dvv;</div>
<div class="line"><a name="l00375"></a><span class="lineno">  375</span>&#160;    cout&lt;&lt;<span class="stringliteral">&quot; merging &quot;</span>&lt;&lt;member.size()&lt;&lt;<span class="stringliteral">&quot; -------&gt; &quot;</span>;</div>
<div class="line"><a name="l00376"></a><span class="lineno">  376</span>&#160;    <span class="keywordflow">while</span>(nid &lt; member.size())</div>
<div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160;    {</div>
<div class="line"><a name="l00378"></a><span class="lineno">  378</span>&#160;        <span class="keywordflow">for</span>(<span class="keyword">auto</span> p = member.begin(); p != member.end(); p ++)</div>
<div class="line"><a name="l00379"></a><span class="lineno">  379</span>&#160;        {</div>
<div class="line"><a name="l00380"></a><span class="lineno">  380</span>&#160;            <span class="keywordflow">if</span>(p-&gt;id == member[nid].id)</div>
<div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;            {</div>
<div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;                <span class="keywordflow">continue</span>;</div>
<div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;            }</div>
<div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;            dvv = member[nid].velocity - p-&gt;velocity;</div>
<div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;            <span class="keywordtype">double</span> delta_energy = abs(member[nid].energy - p-&gt;energy);</div>
<div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;            <span class="keywordtype">double</span> delta_velocity = Mod(dvv);</div>
<div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;            <span class="keywordflow">if</span>(delta_energy / member[nid].energy &lt;= delta_e &amp;&amp; delta_energy / p-&gt;energy &lt;= delta_e &amp;&amp; \</div>
<div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;                    (delta_velocity / Mod(member[nid].velocity) &lt; delta_v &amp;&amp; Mod(dvv) / \</div>
<div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;                        Mod(p-&gt;velocity) &lt; delta_v))</div>
<div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;            {</div>
<div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;                member[nid].energy = (member[nid].energy * member[nid].weight + p-&gt;energy * p-&gt;weight);</div>
<div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160;                member[nid].velocity = (member[nid].velocity * member[nid].weight + p-&gt;velocity * p-&gt;weight);</div>
<div class="line"><a name="l00393"></a><span class="lineno">  393</span>&#160;                member[nid].position = (member[nid].position + p-&gt;position) * 0.5;</div>
<div class="line"><a name="l00394"></a><span class="lineno">  394</span>&#160;                member[nid].weight += p-&gt;weight;</div>
<div class="line"><a name="l00395"></a><span class="lineno">  395</span>&#160;                member[nid].energy /= member[nid].weight;</div>
<div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;                member[nid].velocity = member[nid].velocity / member[nid].weight;</div>
<div class="line"><a name="l00398"></a><span class="lineno">  398</span>&#160;                <span class="comment">//cout&lt;&lt;p-&gt;id&lt;&lt;&quot; will be del and merge with &quot;&lt;&lt;member[nid].id&lt;&lt;endl;</span></div>
<div class="line"><a name="l00399"></a><span class="lineno">  399</span>&#160;                p = member.erase(p);</div>
<div class="line"><a name="l00400"></a><span class="lineno">  400</span>&#160;                p --;</div>
<div class="line"><a name="l00401"></a><span class="lineno">  401</span>&#160;            }</div>
<div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160;        }</div>
<div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;        nid ++;</div>
<div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;    }</div>
<div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;    cout&lt;&lt;<span class="stringliteral">&quot; &quot;</span>&lt;&lt;member.size()&lt;&lt;endl;</div>
<div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;    <a class="code" href="../../d6/d19/class____Ptclslist.html#a682e2a2543c54e627c6fc3fd42843878">squeeze</a>();</div>
<div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;}</div>
<div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;</div>
<div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;</div>
<div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;__Ptcls::__Ptcls(<span class="keyword">const</span> <a class="code" href="../../d0/de9/structParticles.html">Particles</a>&amp; myptcl)</div>
<div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;{</div>
<div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;    <span class="keywordtype">id</span> = <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a>;</div>
<div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">total_id</a> ++;</div>
<div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">mass</a> = myptcl.mass;</div>
<div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">charge</a> = myptcl.charge;</div>
<div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">energy</a> = myptcl.energy;</div>
<div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;    species = myptcl.species;</div>
<div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;    <a class="code" href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">velocity</a> = myptcl.velocity;</div>
<div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;    lorentz_gamma = 1.0 + myptcl.energy / (myptcl.mass * c_light_speed * c_light_speed);</div>
<div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;    </div>
<div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;}</div>
<div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;</div>
<div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="ttc" id="class____Ptclslist_html_a639c350be512a5b80043429e5e873a75"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a639c350be512a5b80043429e5e873a75">__Ptclslist::species</a></div><div class="ttdeci">int species</div><div class="ttdoc">当前粒子种类 </div><div class="ttdef"><b>Definition:</b> particles.h:191</div></div>
<div class="ttc" id="class____Vect2_html"><div class="ttname"><a href="../../d6/d9c/class____Vect2.html">__Vect2&lt; double &gt;</a></div></div>
<div class="ttc" id="class____Ptclslist_html_ab01ad08b2b0eba942ccb4e7af6b05ea6"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#ab01ad08b2b0eba942ccb4e7af6b05ea6">__Ptclslist::get_ptcl</a></div><div class="ttdeci">__Ptcls &amp; get_ptcl(int offset)</div><div class="ttdoc">返回offset处粒子信息 </div><div class="ttdef"><b>Definition:</b> particles.h:259</div></div>
<div class="ttc" id="class____Ptcls_html_a7ef43b6a1b34878da612e73bc643e20d"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a7ef43b6a1b34878da612e73bc643e20d">__Ptcls::weight</a></div><div class="ttdeci">long weight</div><div class="ttdoc">粒子权重 </div><div class="ttdef"><b>Definition:</b> particles.h:76</div></div>
<div class="ttc" id="class____Ptclslist_html_a98e6733c424414f7cc596e64486e43b2"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a98e6733c424414f7cc596e64486e43b2">__Ptclslist::iter</a></div><div class="ttdeci">vector&lt; __Ptcls &gt;::iterator iter</div><div class="ttdoc">粒子类的一个迭代器，一般不要调用，除非你懂 </div><div class="ttdef"><b>Definition:</b> particles.h:189</div></div>
<div class="ttc" id="class____Ptcls_html"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html">__Ptcls</a></div><div class="ttdoc">单个粒子类信息，可以是任何种类粒子 </div><div class="ttdef"><b>Definition:</b> particles.h:9</div></div>
<div class="ttc" id="class____Ptcls_html_af863f11743d992891e0e79de66e28f65"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#af863f11743d992891e0e79de66e28f65">__Ptcls::advance</a></div><div class="ttdeci">void advance(double &amp;delta_t)</div><div class="ttdoc">粒子推动 </div><div class="ttdef"><b>Definition:</b> particles.h:112</div></div>
<div class="ttc" id="class____Ptclslist_html_a7feae0ba57fcfc32edaffa105cb1d09c"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a7feae0ba57fcfc32edaffa105cb1d09c">__Ptclslist::size</a></div><div class="ttdeci">int size()</div><div class="ttdoc">返回当前list中粒子个数 </div><div class="ttdef"><b>Definition:</b> particles.h:247</div></div>
<div class="ttc" id="class____Ptclslist_html_a86188e07e2f77f26cb58aaab4fa92743"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a86188e07e2f77f26cb58aaab4fa92743">__Ptclslist::operator=</a></div><div class="ttdeci">__Ptclslist operator=(const __Ptclslist &amp;other)</div><div class="ttdoc">通过重载可以直接复制了 </div><div class="ttdef"><b>Definition:</b> particles.h:300</div></div>
<div class="ttc" id="class____Ptclslist_html_ac6fc791bb5b3b456b7a909c175f43340"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#ac6fc791bb5b3b456b7a909c175f43340">__Ptclslist::mylength</a></div><div class="ttdeci">__Vect2&lt; double &gt; mylength</div><div class="ttdoc">网格长宽 </div><div class="ttdef"><b>Definition:</b> particles.h:193</div></div>
<div class="ttc" id="class____Ptcls_html_a18c3f527195910edacebf07d1deee7c3"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a18c3f527195910edacebf07d1deee7c3">__Ptcls::lorentz</a></div><div class="ttdeci">void lorentz(double &amp;delta_t)</div><div class="ttdoc">粒子根据速度推动一步 </div><div class="ttdef"><b>Definition:</b> particles.h:117</div></div>
<div class="ttc" id="class____Ptcls_html_a39e0e214e5aaedc38625b1d2d3eab80a"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a39e0e214e5aaedc38625b1d2d3eab80a">__Ptcls::old_pos</a></div><div class="ttdeci">__Vect2&lt; double &gt; old_pos</div><div class="ttdoc">前一步的位置 </div><div class="ttdef"><b>Definition:</b> particles.h:66</div></div>
<div class="ttc" id="class____Ptcls_html_aeba79577ec36efdd5d55f0b483663eeb"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#aeba79577ec36efdd5d55f0b483663eeb">__Ptcls::id</a></div><div class="ttdeci">long id</div><div class="ttdoc">粒子的标号 </div><div class="ttdef"><b>Definition:</b> particles.h:60</div></div>
<div class="ttc" id="class____Ptcls_html_aeff96eb3ac159cdefef7f8e2fd03f49f"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#aeff96eb3ac159cdefef7f8e2fd03f49f">__Ptcls::bfield</a></div><div class="ttdeci">__Vect3&lt; double &gt; bfield</div><div class="ttdoc">粒子所在位置磁场 </div><div class="ttdef"><b>Definition:</b> particles.h:74</div></div>
<div class="ttc" id="class____Ptcls_html_a9df3807cf12885312f760d8c7b400365"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a9df3807cf12885312f760d8c7b400365">__Ptcls::efield</a></div><div class="ttdeci">__Vect3&lt; double &gt; efield</div><div class="ttdoc">粒子所在位置电场 </div><div class="ttdef"><b>Definition:</b> particles.h:72</div></div>
<div class="ttc" id="class____Ptclslist_html_a51c29b9954cc6f3fe9c51fb919d38a6d"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a51c29b9954cc6f3fe9c51fb919d38a6d">__Ptclslist::pop</a></div><div class="ttdeci">void pop(int offset)</div><div class="ttdoc">删除offset处的粒子，一般不会主函数调用，不要使用，除非你懂 </div><div class="ttdef"><b>Definition:</b> particles.h:230</div></div>
<div class="ttc" id="class____Ptclslist_html"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html">__Ptclslist</a></div><div class="ttdoc">一个同类粒子的粒子束，隶属于一个网格，也可以不是 </div><div class="ttdef"><b>Definition:</b> particles.h:182</div></div>
<div class="ttc" id="class____Ptcls_html_ae16bdf32f375a3037b79b977dafb5b29"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#ae16bdf32f375a3037b79b977dafb5b29">__Ptcls::energy</a></div><div class="ttdeci">double energy</div><div class="ttdoc">粒子能量，没有乘以权重 </div><div class="ttdef"><b>Definition:</b> particles.h:78</div></div>
<div class="ttc" id="class____Ptclslist_html_a682e2a2543c54e627c6fc3fd42843878"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a682e2a2543c54e627c6fc3fd42843878">__Ptclslist::squeeze</a></div><div class="ttdeci">void squeeze()</div><div class="ttdoc">压缩空间 </div><div class="ttdef"><b>Definition:</b> particles.h:203</div></div>
<div class="ttc" id="class____Ptcls_html_aec1deca97ac099553420c1f9cc138b3b"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#aec1deca97ac099553420c1f9cc138b3b">__Ptcls::charge</a></div><div class="ttdeci">double charge</div><div class="ttdoc">粒子电荷，没有乘以权重 </div><div class="ttdef"><b>Definition:</b> particles.h:80</div></div>
<div class="ttc" id="class____Ptclslist_html_abd6d139f042a71a595aa7f4aa2a21714"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#abd6d139f042a71a595aa7f4aa2a21714">__Ptclslist::get_sp_ek</a></div><div class="ttdeci">double get_sp_ek(const int &amp;sp)</div><div class="ttdoc">返回该list的能量 </div><div class="ttdef"><b>Definition:</b> particles.h:357</div></div>
<div class="ttc" id="class____Ptcls_html_a4c6caa099428ac9c06c2c96c78456021"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a4c6caa099428ac9c06c2c96c78456021">__Ptcls::old_vel</a></div><div class="ttdeci">__Vect3&lt; double &gt; old_vel</div><div class="ttdoc">velocity previous step </div><div class="ttdef"><b>Definition:</b> particles.h:70</div></div>
<div class="ttc" id="class____Ptclslist_html_aad813503cf966e45b694006e71a022e1"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#aad813503cf966e45b694006e71a022e1">__Ptclslist::showme</a></div><div class="ttdeci">void showme()</div><div class="ttdoc">输出每个粒子 </div><div class="ttdef"><b>Definition:</b> particles.h:294</div></div>
<div class="ttc" id="class____Ptcls_html_acb07ae3c6b8df2cfa5dacdf2aa988ee7"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#acb07ae3c6b8df2cfa5dacdf2aa988ee7">__Ptcls::position</a></div><div class="ttdeci">__Vect2&lt; double &gt; position</div><div class="ttdoc">粒子坐标,全局绝对 </div><div class="ttdef"><b>Definition:</b> particles.h:64</div></div>
<div class="ttc" id="class____Ptclslist_html_a89bd5130d9b94c556e1f23ccf1a98485"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a89bd5130d9b94c556e1f23ccf1a98485">__Ptclslist::clear</a></div><div class="ttdeci">void clear()</div><div class="ttdoc">清除当前list </div><div class="ttdef"><b>Definition:</b> particles.h:241</div></div>
<div class="ttc" id="class____Ptcls_html_a77ff6b30069c0a76772206cebceb009f"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a77ff6b30069c0a76772206cebceb009f">__Ptcls::set_field</a></div><div class="ttdeci">void set_field(const __Vect3&lt; double &gt; &amp;oefield, const __Vect3&lt; double &gt; &amp;obfield)</div><div class="ttdoc">设置粒子感受到的场 </div><div class="ttdef"><b>Definition:</b> particles.h:103</div></div>
<div class="ttc" id="class____Ptclslist_html_abda46652d821fa0b406851572ab9a7a9"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#abda46652d821fa0b406851572ab9a7a9">__Ptclslist::member</a></div><div class="ttdeci">vector&lt; __Ptcls &gt; member</div><div class="ttdoc">粒子类成员 </div><div class="ttdef"><b>Definition:</b> particles.h:185</div></div>
<div class="ttc" id="class____Ptcls_html_a190602549e77f92d18694bfac93774de"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a190602549e77f92d18694bfac93774de">__Ptcls::species</a></div><div class="ttdeci">int species</div><div class="ttdoc">粒子的种类：0电子，1：离子，2：正电子，3：光子 </div><div class="ttdef"><b>Definition:</b> particles.h:62</div></div>
<div class="ttc" id="structParticles_html"><div class="ttname"><a href="../../d0/de9/structParticles.html">Particles</a></div><div class="ttdef"><b>Definition:</b> type.h:373</div></div>
<div class="ttc" id="class____Ptcls_html_ae27fae8bcdf7b0f95d60931480a149d4"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#ae27fae8bcdf7b0f95d60931480a149d4">__Ptcls::total_id</a></div><div class="ttdeci">static long total_id</div><div class="ttdoc">current max id </div><div class="ttdef"><b>Definition:</b> particles.h:56</div></div>
<div class="ttc" id="class____Ptcls_html_a0709048b1598219650298d09ecda7a6b"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a0709048b1598219650298d09ecda7a6b">__Ptcls::velocity</a></div><div class="ttdeci">__Vect3&lt; double &gt; velocity</div><div class="ttdoc">粒子速度 </div><div class="ttdef"><b>Definition:</b> particles.h:68</div></div>
<div class="ttc" id="class____Ptclslist_html_a51c4123e00a360321de6765592ee474e"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a51c4123e00a360321de6765592ee474e">__Ptclslist::weightfunc</a></div><div class="ttdeci">double weightfunc(const __Vect2&lt; double &gt; &amp;mycoords, const __Vect2&lt; double &gt; &amp;other, const __Vect2&lt; double &gt; &amp;cell_length)</div><div class="ttdoc">与cellinfo中相同 </div><div class="ttdef"><b>Definition:</b> particles.h:315</div></div>
<div class="ttc" id="class____Vect3_html"><div class="ttname"><a href="../../da/d1a/class____Vect3.html">__Vect3&lt; double &gt;</a></div></div>
<div class="ttc" id="class____Ptcls_html_ae2c78f356830d91749bfc3cd6f7ffa6c"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#ae2c78f356830d91749bfc3cd6f7ffa6c">__Ptcls::total_num</a></div><div class="ttdeci">static long total_num</div><div class="ttdoc">current box total num </div><div class="ttdef"><b>Definition:</b> particles.h:58</div></div>
<div class="ttc" id="class____Ptcls_html_a4c8ff1d2388163601e00c6685af2859a"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a4c8ff1d2388163601e00c6685af2859a">__Ptcls::mass</a></div><div class="ttdeci">double mass</div><div class="ttdoc">粒子质量，没有乘以权重 </div><div class="ttdef"><b>Definition:</b> particles.h:82</div></div>
<div class="ttc" id="class____Ptcls_html_a522b307a9a3fefda5c1cf489ca77fa9b"><div class="ttname"><a href="../../d9/ddf/class____Ptcls.html#a522b307a9a3fefda5c1cf489ca77fa9b">__Ptcls::boris</a></div><div class="ttdeci">void boris(double &amp;delta_t)</div><div class="ttdoc">boris转动进行电磁加速 </div><div class="ttdef"><b>Definition:</b> particles.h:127</div></div>
<div class="ttc" id="class____Ptclslist_html_abf8f9068e29bee8f2cdef145ba7e32ff"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#abf8f9068e29bee8f2cdef145ba7e32ff">__Ptclslist::advance</a></div><div class="ttdeci">void advance(double delta_t, vector&lt; __Ptcls &gt; &amp;lists, const __Vect2&lt; double &gt; &amp;xymin, const __Vect2&lt; double &gt; &amp;xymax, const int &amp;myrank)</div><div class="ttdoc">推动该list中所有粒子一步 </div><div class="ttdef"><b>Definition:</b> particles.h:272</div></div>
<div class="ttc" id="class____Ptclslist_html_a7969ef614256281a5103a4478e009e46"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a7969ef614256281a5103a4478e009e46">__Ptclslist::push</a></div><div class="ttdeci">void push(const __Ptcls &amp;other)</div><div class="ttdoc">将一个粒子压入当前list </div><div class="ttdef"><b>Definition:</b> particles.h:216</div></div>
<div class="ttc" id="class____Ptclslist_html_a67f83b452b3c0970072c6b8454fe8835"><div class="ttname"><a href="../../d6/d19/class____Ptclslist.html#a67f83b452b3c0970072c6b8454fe8835">__Ptclslist::merge</a></div><div class="ttdeci">void merge(const double &amp;delta_e, const double &amp;delta_v)</div><div class="ttdoc">利用此方法对当前的list进行压缩,如光子是需要压缩的.当光子的能量差异在delta_e之内时会发生压缩.当然可以指定空间距离 </div><div class="ttdef"><b>Definition:</b> particles.h:371</div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated on Tue Jun 16 2015 22:48:37 for hybridMPIC by &#160;<a href="http://www.doxygen.org/index.html">
<img class="footer" src="../../doxygen.png" alt="doxygen"/>
</a> 1.8.9.1
</small></address>
</body>
</html>
